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ABSTRACT 

We report on 450- and 850-//m observations of the interacting galaxy pair, VV114E + W (IC 1623), 
taken with the SCUBA camera on the James Clerk Maxwell Telescope, and near-infrared observations 
taken with UFTI on the UK Infrared Telescope. The system VV114 is in an early stage of a gas- 
rich merger. We detect submillimeter (sub-mm) emission extended over 30" (12 kpc) and find a good 
correlation between the spatial distribution of the sub-mm and CO emission. Both the CO and sub-mm 
emission peak near the reddest region of VV 114 E and extend toward VV114W. The bulk of the sub-mm 
emission resides in the central region showing the largest CO velocity gradients, which is thought to mark 
the kinematic centroid of the merger remnant. We derived a total dust mass of 1.2 x 10 8 M Q , assuming 
a power-law distribution of dust temperatures. The sub-mm observations suggest that the majority of 
the dust is relatively cool (Td ~ 20-25 K), and the total dust mass is about 4 times higher than that 
inferred from the IRAS data alone. The system will likely evolve into a compact starburst similar to 
Arp220. 

Subject headings: galaxies: evolution — galaxies: individual (VV114) — galaxies: interactions — 
galaxies: starburst 



1. INTRODUCTION 

Many extremely luminous infrared starbursts (I/ir > 
10 Lq) are merging systems, containing large reservoirs 
of molecular gas and dust (Sanders & Mirabel 1996, and 
references therein). To understand the evolution of the 
ISM and star formation in these systems, we require imag- 
ing observations of the molecular gas and dust reservoirs 
associated with the star-forming regions. Much of what is 
known about these systems has been inferred from inter- 
ferometric CO observations (e.g., Scoville, Yun, & Bryant 
1997; Downes & Solomon 1998). Our knowledge of the 
distribution of thermal dust emission in merging systems 
is lacking, due to the limited spatial resolution of the In- 
frared Astronomical Satellite (IRAS). With the commis- 
sioning of the new Sub-mm Common-User Bolometer Ar- 
ray (SCUBA, Holland et al. 1999) camera on the James 
Clerk Maxwell Telescope (JCMT), we can now directly im- 
age the thermal dust emission at relatively high (~ 8 — 15") 
spatial resolution. 

The interacting galaxy pair VV 114 is only 80 Mpc away 
(H = 75 km s" 1 Mpc" 1 ) and is comprised of two stellar 
components separated by 15" (6 kpc in projection), desig- 
nated VV 114 E and VV 114 W (Knop et al. 1994). The ex- 
treme infrared colors of VV 114E indicate the presence of 
a large concentration of dust, while VV 114 W is relatively 
unobscured at optical wavelengths (Knop et al. 1994). The 
CO emission is distributed along a bar-like structure lo- 
cated in between the optical components (Yun, Scoville, 
& Knop 1994). The extended nature of the CO emission 
as well as the non-circular gas kinematics suggests that 
VV 114 is in an early stage of a gas-rich merger (Yun et al. 



1994). The gas in VV 114 has already begun to flow into 
its central regions, while the stellar components are still 
well separated. This is consistent with theoretical models 
of young mergers (Barnes & Hernquist 1991, 1996). 

In addition to the capability of SCUBA to resolve 
the VV114 system, the sub-mm observations at 450 and 
850 /im allow us to constrain the amount of cool mate- 
rial present (Td <; 30 k). Preliminary evidence for the 
existence of cold dust in galaxies comes from the gas- 
to-dust ratios estimated using the warm dust masses de- 
rived from IRAS. The typical gas-to-dust ratios found for 
the bright IRAS galaxies are about an order of magni- 
tude larger than the Galactic value derived from extinction 
studies, which includes dust at all temperatures (Devereux 
& Young 1990). A simple explanation for this discrepancy 
is that the majority of the dust in bright IRAS galaxies 
is cooler (Td <J 30 k) than the warm material observed by 
IRAS. Observations with SCUBA can directly search for 
this component and provide a test of this hypothesis. 

2. OBSERVATIONS 

2.1. Sub-mm Imaging 

Fully-sampled images of VV 1 14 were obtained using 
the SCUBA sub-mm camera during 1997 July and 1998 
May. We used the 91-element short-wave array at 450 /im 
(7"8 FWHM beam) and the 37-element long-wave array at 
850 /im (14'.' 7 FWHM beam) during excellent conditions. 
Typical 850-/im zenith opacities were approximately 0.2, 
determined from hourly skydips. We obtained 7.7 ks of 
useful on-source integration time at 850 /im and 5.1 ks at 
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Fluxes were calibrated against Uranus and the flux den- 
sities determined for VV114 at 450 and 850 /im are ac- 
curate to ±18% and ±14%, respectively, including errors 
induced by opacity uncertainties and the absolute flux un- 
certainty of Uranus. The uncertainty in the 450/im/850/im 
flux ratio is approximately ±7%, which represents the 
full dispersion in the measurements made on the separate 
nights. The error in the 450^m/850^m flux ratio is less 
than total uncertainty at either wavelength since the abso- 
lute flux calibration error of Uranus cancels for the ratio. 
Pointing was checked regularly using a nearby blazar. We 
estimate the overall positional accuracy of the submillime- 
ter maps to be ±2". 

2.2. Near-infrared Imaging 

Near-infrared, J- and if-band, imaging of VV114 was 
obtained during the night of 1998 October 19 using the 
new UKIRT Fast-Track Imager (UFTI) on the 3.8-m UK 
Infrared Telescope (UKIRT) on Mauna Kea. UFTI com- 
prises a 1024 2 HgCdTe HAWAII array cooled to 77 k and 
sensitive from 0.8 to 2.5/zm. The pixel scale is 0'.'091, al- 
lowing Nyquist sampling of the best seeing experienced on 
UKIRT, and the field of view is 92". The observations of 
VV114 consisted of three sets of 9 exposures, two in K 
(2.2/um) and one in J (1.2/Lim), each of 30s, dithered on a 
3x3 grid with 9'/l spacing. These were interspersed by 
similar length observations of an offset sky region and dark 
exposures to allow the removal of the sky and instrumental 
backgrounds. The final stacked exposures represents 270 s 
in J and 540 s in K. The tip/tilt secondary on UKIRT 
was used to provide standard adaptive correction within 
our field, although the conditions were non-photometric 
and the seeing achieved was only mediocre, ~ Of! 7. Never- 
theless, these images represent a substantial improvement 
in resolution over those presented by Knop et al. (1994). 

3. RESULTS 

Figures la & lb show the 450- and 850-/zm images of 
VV 114 taken with SCUBA. The integrated flux densities 
summed over the emission regions are 2.43 ± 0.44 Jy and 
0.273 ± 0.038 Jy at 450^m and 850^m, respectively (Ta- 
ble 1). At 850/xm, we detect emission extended in the east- 
west direction, while at 450^m we appear to resolve the 
central emission into two peaks (450/zm[l] and 450^m[2] 
in Table 1). Given that the separation of the two 450/im 
peaks is near the resolving limit of SCUBA, observations 
at higher resolution are required to confirm this double- 
peaked morphology. When the 450-/xm data are convolved 
to the resolution of the 850-yum map, the data at both 
wavelengths show a similar shape, size, and position. At 
this resolution (14'.'7), we fail to detect any significant vari- 
ations in the S'(450^m)/S , (850^m) flux ratio. Integrating 
over all emission regions, we find a flux ratio of 8.9 ± 0.6, 
which is also the ratio measured for the peak of the sub- 
mm emission (at a resolution of 14'.'7). The consistencies in 
the 450-/im and 850-^im data would be expected if the bulk 
of the sub-mm emission is on the black-body tail of spec- 
tral energy distribution with dust temperatures T ;> 10 K. 
Observations at higher resolution are required to search 
for possible variations in the S'(450/im)/S , (850/im) ratio 
on smaller spatial scales. 

The sub-mm emission lies along the CO bar-like struc- 
ture (Fig. lc), extending east- west between the two opti- 



cal components. Both the CO and sub-mm emission peak 
near the brightest K— band source of VV 114 which is by 
far the reddest part of the VV 114 system, as indicated by 
the J - K color map (Fig. Id). Although the NIR data 
were not taken in photometric conditions, an approximate 
J — K magnitude scale was derived by matching the data 
with the previously determined J — K color of VV 114W 
(Knop et al. 1994). The higher-resolution data presented 
here suggests a color of ( J — K) ~ 3 for the bright compact 
red component in VV114E. These results are consistent 
with those found for the Near-Infrared Camera and Multi- 
Object Spectrometer (NICMOS) data taken with the Hub- 
ble Space Telescope (Scoville et al. 1999). 

4. DISCUSSION 

4.1. Comparison of the Sub-mm Maps with Data at 
Other Wavelengths 

The sub-mm, radio continuum, and CO emission from 
VV 114 all have roughly the same spatial extent. However, 
there are significant differences in their detailed structure. 
In particular, the radio sources (Condon et al. 1990; Con- 
don et al. 1991) are located near the near-infrared (NIR) 
emission peaks (Knop et al. 1994), suggesting that these 
regions are responsible for the majority of the ongoing star 
formation. The sub-mm and CO peaks do not appear to 
be correlated with the radio and NIR peaks but tend to 
lie primary in between VV 114E + W. 

The CO peak lies approximately 2" from the submil- 
limcter peak, which is consistent within the positional un- 
certainties of the data sets. The general similarities of the 
CO and sub-mm maps suggest that both the CO(l— >0) 
and sub-mm emission probe the same material: the molec- 
ular gas and dust reservoir associated with the merger 
event. Despite these consistencies, there are differences in 
their detailed structure which may indicate optical depth 
or excitation effects in the CO emission. The CO emis- 
sion is distributed more smoothly along the bar-like struc- 
ture between VV114E and VV114W, while the sub-mm 
emission is more strongly peaked near VV 114E. The tidal 
tails are also more apparent in the CO map, while the 
450/im emission is more extended in the east and south- 
west directions. The excess sub-mm emission in the east 
and south-west directions may be associated with Hi tidal 
debris infalling into the central regions (Hibbard & Yun 
1999). 

VV 114E is located near a peak of the sub-mm emission 
which is consistent with the high level of dust obscuration 
inferred from its extremely red NIR colors (Fig. Id; also 
see Knop et al. 1994). VV114W, on the other hand, lies 
along a line of sight which is nearly free of sub-mm emis- 
sion (Fig. la) and is consistent with the low level of ex- 
tinction estimated at NIR wavelengths (Knop et al. 1994). 
The NIR data show that the star formation is occurring 
globally throughout the system, but the most active site 
is in the nucleus of VV 114E, which is associated with the 
brightest radio sources in the 8.4-GHz sub-arcsec resolu- 
tion map (Condon et al. 1991). The ISO 15-yum continuum 
also peaks strongly on the eastern nucleus of VV 114 (Yun 
et al. 1999). VV114E, itself, is comprised of two bright 
nuclear sources in the NIR (Knop et al. 1994). The north- 
east component of VV 114 E is strongest in J and is asso- 
ciated with the brightest compact radio source (marked by 
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a triangle symbol in Fig. 1) in the high-resolution radio 
map (Condon et al. 1991). The red south-west component 
of VV114E is strongest in K (marked by a x symbol in 
Fig. 1) and is associated with a region containing sev- 
eral bright, compact radio sources. This red south-west 
component of VV114E is nearest to the position of the 
peaks in the sub-mm and CO emission, but is still offset 
by 3'.'5-5". This positional offset is larger than the posi- 
tional errors in the data sets and hence may be significant. 

Although most of the current star formation is occur- 
ring in and around the stellar components seen in the NIR 
(Knop et al. 1994), the central sub-mm and CO emission 
regions between the optical components may mark the lo- 
cation of a future major starburst given that these regions 
show the largest CO velocity gradients (Yun et al. 1994) 
and contain the bulk of the molecular gas and dust in 
the system. The apparent displacement of the sub-mm 
and CO emission from the most active regions of star- 
formation is consistent with the low HCN/CO ratio ob- 
served in VV114 (Gao 1996). Typical infrared luminous 
galaxies show enhanced HCN/CO ratios, presumably due 
to the high density of the molecular gas in the starburst re- 
gions (Gao 1996). The low HCN/CO ratio in VV 114 may 
indicate that the bulk of the molecular gas is in regions of 
lower density and is displaced from the most active star- 
forming regions. Although the sub-mm emission of VV 114 
is currently extended over 12 kpc, it is likely to evolve into 
a more compact starburst system similar to that seen for 
Arp 220 (Downes & Solomon 1998; Sakamoto et al. 1998). 
In the context of the numerical models for merging galax- 
ies (Barnes & Hernquist 1996), VV114 is an early-stage 
merger where the gas is accumulating within the central 
regions in advance of the stars. The system has yet to 
trigger a strong compact starburst at the dynamical cen- 
troid as seen for Arp 220 and other compact ultraluminous 
galaxies (Downes & Solomon 1998). 

4.2. Estimating the Total Dust Mass 

Observations at submillimeter wavelengths provide a 
better estimate of the total dust mass than that inferred 
from the IRAS data alone, since IRAS was not sensi- 
tive to cool dust. Assuming optically thin dust emission, 
the mass in dust is given by M d = S v D 2 (fed B[v, Td]) -1 , 
where B(v, Td) is the black-body function for a frequency 
v and a temperature Td, D is the distance, and fed is 
the dust absorption coefficient. We adopt fed = 10 g _1 
cm 2 (A/250/rni)"' 3 , where (3 = 1 for A < 25(tyan (Hilde- 
brand 1983). At longer wavelengths (A > 250^m), the 
typical values for fi range between 1 and 2. Empirically 
at long wavelengths, the Galactic cirrus is well fitted by 
(3 ~ 1.5 (Masi et al. 1995), while Galactic star-forming 
regions have (3 ~ 2 (Chini et al. 1986; Lis et al. 1998). 

From the submillimeter measurements, we can directly 
estimate the value of j3 appropriate for VV 114 as a func- 
tion of temperature. The integrated sub-mm flux density 
ratio is S(450/zm)/S(850/Ltm) = 8.9 ± 0.6. The value for (3 
is calculated from 5'(450^m)/S'(850Mm) = (450/850)-^ x 
5(666 GHz, T d )/S(353 GHz, T d ). At T d = 20k, the data 
imply (3 = 2.15 ± 0.11. At the temperature of T d = 41 
inferred from the IRAS data (Soifer et al. 1989), the sub- 
mm data imply (3 = 1.76 ±0.11. Since the mass-weighted 
average temperature of the dust is expected to be less than 
or equal to the IRAS dust temperature, the lower limit for 



(3 in VV 114 is (3 > 1.65. In fact, we find (see below) that 
most of the dust is relatively cool in VV 114 (T ~ 20-25 k) 
suggesting that (3 ~ 2, which is similar to that observed 
for Galactic star-forming regions (Chini et al. 1986). 

We fit thermal dust spectral energy distributions (SEDs) 
to the observational data in order to estimate the total 
dust mass in VV 114. Figure 2a shows models appropriate 
for a temperature of Td — 41 ± 4 K derived from the IRAS 
60^m/100^m flux density ratio (Soifer et al. 1989). The 

dust mass inferred for this temperature is (3± 2 ) x 10 7 M Q 
(Table 2). The uncertainty in fed provides an additional 
uncertainty of about a factor of two in the dust mass. 
None of the single-component temperature models fit the 
sub-mm data. The low (3 models are inconsistent with the 
5'(450^m)/S'(850^m) flux density ratio, and for (3 ~ 2 the 
sub-mm fluxes are significantly higher than those expected 
from the IRAS dsXa,. A simple interpretation of the data is 
the presence of cooler material which was not detected by 
IRAS. Figure 2b shows a three-component fit to the data, 
assuming (3 — 2 at A > 250^m. The mass of dust residing 
in the cool component (T d = 20k) is 1-2 x 10 8 M Q , which 
is significantly larger than that found in the warm IRAS 
component (T d = 41 K, 0.3 x 10 8 M ). By varying the al- 
lowed dust temperatures, the relative amount of dust in 
the warm (37-45 k) and cool (~15-25 k) components can 
be modified by factors of 2-3. The amount of dust in the 
hot component (Td = 150 k), constrained by the 12 and 
25^m data, is negligible (Table 2). 

More complicated models, which include several differ- 
ent grain types, cool cirrus, and warm starburst compo- 
nents could also be fitted to the data (Rowan-Robinson 
1992; Andrcani & Franccschini 1996). Although such mod- 
els may be more realistic, the number of model parame- 
ters exceeds the number of observational data points for 
VV114. For the model with three temperature compo- 
nents (Fig. 2b), there are six free parameters (mass and 
temperature of each component) for only six data points. 
In reality, there is likely to be a continuous range of dust 
temperatures. 

A simple approach is to adopt a power-law distribution 
of dust temperatures; dM^/dT oc T~ x (Xie, Goldsmith, 
& Zhou 1991). The total dust mass is determined by in- 
tegrating the distribution function over a range of tem- 
peratures from T ow to Tiigh- The upper limit to the dust 
temperature has little effect on the total dust mass (e.g., de 
Muizon & Rouan 1985). We assume T h i g h = 200 K. Fig. 2c 
shows dust SEDs computed for several different tempera- 
ture distributions. Considering the small number of free 
parameters (x, T\ OWl and M d [total]), the shape of the com- 
puted SEDs match the observational data remarkably well. 
The best fit to the 60-850^m observational data occur for 
Ti ow = 22 and x ~ 6. The 12-25 /im data are not fitted 
well with a single power-law distribution, but these data 
have little effect on the total dust mass. Since Md oc T^ 5 
for dust in thermal equilibrium (Soifer et al. 1989), we ex- 
pect x ~ 6. With this piece of theoretical insight, the SED 
is effectively fitted by only T\ ow and M d (total), after as- 
suming an appropriate form for £; d (A). The best fit (x = 6; 
T low = 22 k) implies a total dust mass of 1.2 x 10 8 M Q . As- 
suming a power-law distribution of dust temperatures, the 
total dust mass is well constrained by the data. The er- 
ror in fitting the dust mass is only about 25%. However, 
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considering the uncertainty in k& at sub-mm wavelengths 
(Hughes, Dunlop, & Rawlings 1997), the derived dust mass 
is uncertain by approximately a factor of two. 

4.3. Implications for Cool Dust in Luminous IRAS 
Galaxies 

Most of the dust in VV114 is at a temperature of 
T d =20-25 K. Based on the SCUBA data, the total dust 
mass derived for VV114 is four times larger than that 
derived from the 60-100^m IRAS data. Since the IRAS 
colors for VV114 are typical for its infrared luminosity 
(Soifer & Neugebauer 1991), we could expect to find sim- 
ilar results for other luminous IRAS sources. The exis- 
tence of cool dust (Td < 30 k) undetected by IRAS is not 
a surprising result. As stated in §1, the high gas-to-dust 
ratios in galaxies derived from the IRAS data of spiral 
galaxies imply the presence of cool material (Devereux & 
Young 1990). Millimeter observations of spiral galaxies 
also suggest massive cool dust reservoirs with Td ~ 10- 
20 k (Guelin et al. 1993, 1995; Franceschini & Andreani 
1995). Recent SCUBA observations of NGC891 support 
these results with the detection of a cool dust reservoir 
which is over an order of magnitude more massive than its 
warm dust component (Alton et al. 1998). 

The results for VV114 indicate that even very lumi- 
nous IRAS galaxies can have cool dust reservoirs. Based 
solely on the IRAS data, the gas-to-dust ratio for VV 114 
is M(H 2 )/M d = 1100 (Sanders, Scoville, & Soifer 1991), 
which is much higher than the Galactic value of about 100 
(Devereux & Young 1990). By using the dust mass im- 
plied by the sub-mm data, we find a more realistic value 
of Al (H2)/Afd ~ 300. These gas-to-dust ratios assume the 
Galactic CO to H2 conversion factor. If the Galactic value 
is not applicable for ultraluminous galaxies (Solomon et 
al. 1997), then the gas-to-dust ratio for VV114 could be 
reduced further by a factor of 2-4. Observations of other 
luminous IRAS galaxies will test the generality of these 
results. 

The cool dust temperature for VV 114 may be related to 
its early-merger evolutionary phase. More evolved, lumi- 
nous mergers have warm dust emission associated with a 
compact starburst and/or an active galactic nucleus (Maz- 
zarella, Bothun, & Boroson 1991). For example, the sub- 
mm data for the evolved, compact starburst, Arp 220 



(Rigopoulou, Lawrence, & Rowan- Robinson 1996), show 
no evidence for excess sub-mm emission associated with 
cool dust and can be fitted with a single-temperature, 
warm dust model of T d = 47k (Klaas et al. 1997). The 
dust may be cooler in VV 114 because the majority of the 
dust is displaced from the star-forming regions traced by 
the NIR and radio emission peaks. Given the radio-to- 
FIR relationship (Helou, Soifer, & Rowan-Robinson 1985), 
we expect to find warm dust within the star-forming re- 
gions, especially near the most active region, VV114E. 
As VV114 evolves into a compact gaseous system similar 
to Arp 220, the star formation is expected to increase in 
the central regions and the dust temperature should rise 
accordingly. 

5. CONCLUSIONS 

We present sub-mm maps of the young merger system 
VV 114. We detect sub-mm emission in excess of that ex- 
pected from the IRAS data. This sub-mm excess suggests 
the presence of a cool, massive component of dust. By 
fitting a variety of dust models to the SED of VV 114, we 
derive a total dust mass of approximately 1 x 10 s M©, a 
temperature of Td ~ 20-25 K, and (3 ~ 2. The majority 
of the dust is located in between the optical components, 
VV 114E + W, near the suspected dynamical center of the 
merger remnant. The sub-mm emission regions correlate 
well with the CO emission regions, but do not correlate 
well with the peaks in the radio and NIR emission regions 
which are thought to trace the star-formation activity. The 
fact that the sub-mm emission is displaced from the most 
active regions of star formation may explain the cool tem- 
perature for the majority of the dust. Given the extremely 
large reservoir of gas and dust available in the central re- 
gions, VV 114 is expected to evolve into a more luminous 
central starburst, similar to that seen for Arp 220. 
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Fig. 1. — SCUBA (a) 450-^im and (b) 850-/im maps. The Icr rms levels are 75 and 4.4 mjybeam for the 450- and 850^tm data, 
respectively. Contour levels are plotted at a X (-2, 2, 3, 4, 5, 6, 8, 10, 12) for the 450-£tm map and a X (3, 4, 5, 6, 8, 10, 12, 15, 20, 25, 30, 35) 
for the 850-/im map. The 2^tm NIR peaks, signifying the positions of VV114E + W are shown by the X and + symbols (Knop ct al. 1994). 
The 8.4-GHz radio peak is shown by the A symbol (Condon et al. 1991) and the O symbols show peaks in the optical ij-band emission (Knop 
et al. 1994). Panel (c) shows the integrated CO(l— >0) emission (Yun et al. 1994). Contour levels are plotted at 1 Jykm s _1 beam - 1 X (5, 
10, 15, 30, 50, 70). For panels (a), (b), and (c) the beam sizes are shown in the lower left of each panel. Panel (d) shows the NIR (J — K) 
color map, where lighter pixels indicate redder colors and darker pixels represent bluer colors. The color map ranges from (J — K) = 0.9 for 
the blue regions of VV 114 W to the very red compact region in VV 114 E which has a color of (J — K) = 3.0 at the observed resolution (0'.'7). 

Fig. 2. — Model fits to the spectral energy distribution of VV114. The IRAS data points are from Soifer et al. (1989). Panel (a) shows 
single temperature models (Tj = 41k) based on the IRAS 60-100 fim data for f3 = 1, 1.5 (dotted lines), and (3 = 2 (solid line). No single 
temperature model can fit the data. Panel (b) shows a three-component model fit to the data, assuming f3 = 2 for A > 250/im. The three 
different components are plotted as dotted lines, and their sum is shown by the solid line. Panel (c) shows models based on a distribution 
of dust temperatures; dM^/dT oc T~ x (see §4.2). The solid line represents the best model fit with Ti ow = 22 K and x = 6. The dotted lines 
show models with different values of x, while the dashed lines show fits based on different values of Ti ow . 
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Table 1 

submillimeter observational results 



Parameter 


450 y um[l] 450^m[2] 


850/zm 


a(J2000) a ' b 01 h 07 m 


47!55 47!15 


47!35 


5(J2000) a - b -17°30' 


28'.'9 26'.' 1 


27'.' 2 


Beam (FWHM) 


7'.'8 


14'.'7 


Flux (Jy) c 


2.43 ± 0.44 


0.273 ±0.038 



a Positional uncertainty is ±2'.'0. 
b Dcrivcd from a Gaussian fit to the peak. 
c Total integrated flux. 



Table 2 
Dust Models 



Model 


Fitted Data 


?d(K) 


M d (10 7 M Q ) 




Notes 




Single Component 


60-100^m 


41 


3.0 




Fig. 2a 




Three Component 


12-850/um 


20 


13 




Fig. 2b 








41 


2.7 




Fig. 2b 








150 


0.0018 




Fig. 2b 




Temperature Distribution 


450-850^m 


15 


21 


T d 


= T low , x = 6, Fig. 


2c 




450-850^m 


30 


8 




= Ti ow , x = 6, Fig. 


2c 




60-850/im 


22 


12 


T d 


= Ti ow , x = 6, Fig. 


2c 



Frayer et al. 
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Fig. 1 — 
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